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Nontopological solitons are stable field configurations which may be formed in a primor- 
dial phase transition. We study their cosmic evolution and examine the possibility that 
such objects could contribute significantly to the energy density of the Universe. As the 
Universe cools, initially all but the largest lumps evaporate into free particles; those which 
survive may subsequently enter a brief accretion phase before they ‘freeze out’ at a final 
size. Although the minimum critical charges which survive depend on particle masses and 
couplings, we develop an analysis which applies to a wide class of models. In most cases, 
solitons of moderate size survive the evaporation process only if there is a significant charge 
asymmetry or if they form at a temperature well below their binding energy per charge. 
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I. Introduction 

Non-topological solitons (hereafter, NTS’s) are solutions of classical field theories 
which are stable by virtue of a conserved Noether charge carried by fields confined to 
a finite region of space. The prototype for such structures is the phenomenological bag 
model for hadrons, in which light quarks are confined in a false vacuum region by Yukawa 
coupling to a scalar field . 1 In recent years, the bag model has genetically mutated into 
a rich variety of species, which include Q balls , 2 abnormal nuclei , 3 strange quark matter 
nuggets , 4 cosmic neutrino balls , 5 as well as scalar 8 ’ 7 ’ 8 and gauge field 9 non-topological soli- 
tons. Despite the proliferation of models, they all share a number of family resemblances 
which allow them to be classified and discussed in general terms. The key feature of all 
models is that, for particle number Q greater than some minimum Q m in, the confined 
soliton state has lower energy than the free particle state. 

Recently, it has been suggested that non-topological solitons may be abundantly pro- 
duced in a phase transition in the early universe . 4 ’ 7 ' 8 ’ 10 However, as pointed out by Alcock 
and Farhi 11 and Alcock and Olinto 12 in the context of strange matter, at sufficiently high 
temperature, such structures are in general unstable to evaporation into free (unconfined) 
particles. Non-topological solitons are states of low energy but high order. At tempera- 
tures T > Iq, the NTS binding energy per particle, free particles have lower Helmholtz 
free energy, F = E — TS, than the NTS state, due to the entropy term. (This holds as long 
as the temperature is not very near the critical temperature for the phase transition. 10 ) 
At lower temperatures, the energy term dominates, and NTS’s are the preferred states. 

If we focus on the evolution of a single soliton of initial charge (or particle number) 
Q, as the universe cools the lump begins to evaporate. If it is large enough to survive 
down to temperatures T < Iq, it may then enter an accretion phase. Eventually, the 
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accretion or evaporation rate becomes negligible compared to the expansion rate, and 
the system ‘freezes out’ with a fixed charge. Thus, all solitons with initial charge less 
than a critical number Q , die before they reach the accretion or freeze-out phase. Those 
with charges Q Q, are born frozen, i.e., their evaporation and accretion rates are 
always negligible. However, evaporation is so efficient that Q, is often larger than the 
chaxge inside the particle horizon at the time of formation; in these cases, it is unlikely 
that any solitons formed in the phase transition survive to the present. At the other 
extreme, if NTS’s form at a temperature well below Iq, then we expect Q t ~ Qmin • In 
addition, it has recently been suggested that, 13 if the particles are sufficiently strongly 
interacting to maintain chemical equilibrium at temperatures significantly below Iq, then 
a new population of solitons may be built up via accretion and fusion processes at low 
temperature. However, this latter possibility is only feasible if Qmin is very small, of 
order unity (and only if the mass scale is less than about 10 4 GeV). Otherwise formation 
of new solitons by free particle collisions is statistically suppressed at low temperature, 
and the system cannot reach complete chemical equilibrium. Since, in almost all models 
considered in the literature, Q m i n 3> 1, we have relegated discussion of this possibility to 
an Appendix. 

In this paper, we follow the cosmological evolution of non-topological solitons from 
birth, through the adolescent struggle for survival, to placid old age. In the next section, 
we outline the basic properties of NTS’s which we will need; a wide variety of models 
are encapsulated in a small number of parameters. For completeness, in Section III, we 
review several of the NTS formation mechanisms which have been proposed. In Section 
IV, we set up the formalism for treating NTS evolution. We discuss the evolution of single 
NTS lumps in Section V; this treatment is relevant if Q m in is not close to one, in which 
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case the system generally does not reach chemical equilibrium. In Section VI, we discuss 
NTS destruction by ‘boiling’, i.e., bubble nucleation, rather than surface evaporation. Our 
conclusions follow in Section VIII. 

Throughout most of the paper, we shall assume the universe carries a net charge 
asymmetry, and that the particles involved are long-lived. In the subsequent paper of 
this series (hereafter, paper II, in preparation), we drop the assumption of asymmetry 
and discuss how NTS’s evolve as random charge fluctuations in a universe with zero net 
charge. There, we also consider models with unstable particles and discuss cosmological 
constraints on particle physics models with NTS’s. 

II. Properties of Non-topological Solitons 

In all NTS models, a partially ‘confined’ bose or fermi field <f>, carrying a conserved 
additive quantum number, couples to a scalar field tr; the vacuum expectation value o- 0 of 
the ‘confining’ field generates part of the <f> mass, e.g., m $ = m 0 +m^((T 0 ). (I n the case of 
Q-balls, 2 a single complex scalar field performs the roles of both <f> and a.) For appropriate 
couplings and sufficiently large charge, Q > Qmini the lowest energy configuration consists 
not of free <f> particles roaming the vacuum, but of a spherical ‘false vacuum’ region in which 
<r ~ surrounded by a domain wall where a rapidly approaches its ground 6tate value. In 
the NTS interior, < m^(<ro), and the effectively light <f> particles are trapped by the 

mass gap at the domain wall. Since, classically, a non-zero charge implies a time-dependent 
field, the confined particles carry kinetic energy, which varies as an inverse power of the 
volume. The potential (volume) and gradient (surface) energies vary with a positive power 
of the volume. As a consequence, for fixed charge Q, the total energy is minimized at 
a finite radius R , and the resulting configuration is in hydrostatic equilibrium. We will 
consider solutions which carry no currents, in which case the lowest energy configurations 
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axe spherically symmetric. 

For large charges Q > Qmin, the masses of all NTS models follow a simple scaling 

law: 

M(Q) = am+Q*. (1) 

For NTS’s with non-zero volume energy (bag constant), the radius scales as 

R(Q)=(3m^Q^ 3 (volume). (2) 

On the other hand, models in which the interior potential energy density is degenerate 
with the vacuum, U(cr') = U(a o), are confined by surface tension alone; 5,14 in this case, 

R(Q) = /3 , m'^ 1 Q p / 2 (surface). (3) 

In this paper, we will focus almost exclusively on models with non-zero volume energy, 
t.e., we will assume the radius obeys Eqn. (2). The dimensionless parameters p, a, and /? 
depend on coupling constants in the theory. We shall only assume that they satisfy the 
criteria of NTS stability. The stability conditions are 


M(Q) < m$Q 


( 4 ) 


and 


cPM 

dQ 2 


< 0 . 


( 5 ) 


Eqn. (4) ensures stability against decay into free <f> particles, while Eqn. (5) expresses 
stability against fission into smaller soliton fragments. Eqns. (1), (4), and (5) imply the 
upper bound p < 1, and that the minimum stable charge is approximately Qmin — 

(for p < 1). The binding energy per unit charge is defined as the energy required to remove 
a <f> particle from a soliton, 


I Q =m+ + M(Q) - M(Q + 1) 


( 6 ) 
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Using Eqn. (1), we can express this in terms of the mass of the free <f> particles as 

/(<?) = I q /m t = 1 - a((<? + 1)" - Q”] = 1 - api (7) 

Note that for p not very close to one, f(Q) ~ 1 — p(Q/£?mm) p-1 > so f(Q) -»1 as Q/Q m in 
grows large. 

Two notes are in order about the scaling laws, Eqns. (1) and (2). First, although 
they are strictly valid only in the limit Q » Qmin » we will use them as approximate 
formulae for all Q larger than Q m ,„. Second, these scaling laws break down at very large 
charge, Q m ax ~ [(/3 /o:)(m.j>( /^rr^) 2 ] 3 / 2 ^, when gravity becomes important. Since we are 
interested in only moderately sized lumps, Q <C Qmax, we can ignore this complication; 
for discussions of soliton stars, we refer the reader to ref. 15. 

Although the analysis we develop holds for general p, we will often focus on two 
particular cases, p = 1 and 3/4, to demonstrate how the results vary with p. In the 
literature, many models have p essentially equal to one, e.g., Q-balls, fermion NTS’s, 
interacting scalar NTS’s, bag models, and strange matter. In this case, the requirement 
of stability is just a < 1, and the binding energy per charge is independent of Q, i.e., 
f(Q) = 1 — a. (Actually, in these models, due to the surface energy correction, p is slightly 
smaller than one; this ensures stability against fission.) For example, for strange matter, 
f(Q) ~ 0.05 and = m n ~ 1 GeV. Configurations with p = 3/4 arise in a scalar 
model in which the confined field has little or no self-interaction. 6 ’ 7 ’ 8 We shall see that the 
qualitative evolution of NTS’s in the early universe is sensitive to both p and Qmin » as 
well as the mass scale m#. 
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III. Primordial Formation of NTSs 


One can envision several mechanisms by which non-topological solitons may be formed 
in the early universe. In principle, such objects arise in a symmetry-breaking phase tran- 
sition if suitably large regions carrying a net charge Q > Q m in become trapped in the 
false vacuum state. In the simplest scenario, the NTS interior <r = cr' corresponds to a 
metastable local minimum of the potential U(cr) (dashed curve, Fig. 1), and (i) may or (ii) 
may not be a phase of restored symmetry. In case (i), the universe at high temperature 
is initially in the ‘NTS vacuum’ and undergoes a first order phase transition to the true 
vacuum at a critical temperature T c ~ <r 0 . Charge prefers to live in the false vacuum, 
where the <j> mass is smallest. Thus, as bubbles of the true vacuum phase nucleate, grow 
and eventually percolate, the charge may become highly concentrated in the shrinking 
‘NTS vacuum’ regions. Such a mechanism was proposed by Witten for the formation of 
strange matter nuggets in the QCD phase transition. 4 

If the false vacuum is not a phase of restored symmetry, t.e., case (ii) above, then NTS’s 
may form in a second order phase transition. 7 In this instance, during the transition, a 
fraction of space evolves to the metastable minimum <r' instead of the true vacuum (To . 
The relative probability of ending up in the NTS phase is given roughly by the Boltzmann 
formula, PfaUe/ptrue = exp(—8^A f/Ta)- Here, To is the Ginzburg temperature, the 
temperature at which thermal fluctuations between the true and false vacua freeze out, 
(a is the correlation length at the Ginzburg temperature, and Af is the difference in free 
energy density between the true and false vacua. We may interpret Tq as the temperature 
at which NTS’s are first formed; typically To is just below T c , the critical temperature 
for the phase transition. The size distribution of the proto-NTS false vacuum regions is 
calculable from percolation theory, while the scale size is set by the correlation length 
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£g. Generally, false vacuum regions with sizes R £g are exponentially suppressed; 
NTS’s with very large charges are correspondingly rare. This holds particularly if the 
probability pf a he £■ 1/ 3, when only finite NTS regions form. On the other hand, if 
Pfaite > 1/3, both the true vacuum and the NTS phase initially percolate, forming a 
web of infinite interconnecting segments. In this percolating-NTS case, the infinite false 
vacuum regions eventually pinch off into finite blobs when the pressure squeezing them 
from the true vacuum becomes dynamically important. If the vacuum asymmetry is small, 
i.e., A / <£ Tq, pinch-off may occur well after the phase transition is completed; the finite 
NTS’s formed in this way are thus born at a temperature T p i nc h 2g, and have a greater 
chance of surviving. 

In a large class of models, the false vacuum NTS interior is not a metastable minimum 
of the potential, but is instead an unstable, symmetry preserving local maximum (solid 
curve, Fig. 1). In this case, solitons cannot form simply by charge piling up in a metastable 
phase. To see how formation occurs in this case, 10 consider the structure of an NTS (of fixed 
charge) as a function of increasing temperature. Due to finite temperature quantum effects, 
when the temperature reaches the range near Tq, the NTS radius expands enormously. 
If we now imagine starting with a dilute gas of solitons at T ~ 0 and similarly dialing 
up the temperature, then at T ~ Tq, the gas percolates. Reversing the film, as the 
temperature drops below Tq , the Universe, initially in the ‘NTS vacuum’, shatters into 
soliton fragments, which subsequently shrink as they cool. (We note that this process may 
happen at a modest temperature compared to <7o: a large charge asymmetry can delay the 
transition, since charge prefers the NTS vacuum.) The essential idea is that, very close to 
Tq, the free energy of the NTS state becomes small compared to the state of a gas of free 

As a result, solitons can be formed at little cost. Thus, at temperatures T ~ Tg, we 
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expect a population of NTS’s in thermal equilibrium, with an exponential size distribution. 
(Although the initial NTS geometries are expected to be non-spherical, relaxation due to 
surface tension occurs rapidly.) 

From this discussion, several general features emerge. In the scenarios outlined above, 
charge concentration into solitons may arise either due to a net cosmological charge asym- 
metry (like the baryon asymmetry), or due to random fluctuations in a Universe containing 
zero net charge. We will treat the first case in this paper. In addition, if thermal equilib- 
rium is approximately maintained down to a temperature of order Tg, as above, we expect 
a quasi-exponential size distribution of solitons produced in a phase transition. 

Finally, independent of what may happen in a phase transition, NTS’s may form at 
temperatures below Tq by several mechanisms. First, as mentioned in the Introduction, in 
models with small Qmim solitons may form via accretion and fusion after the temperature 
falls below the binding energy Iq. Second, some field theory models give rise to both 
topological and non-topological solitons. For example, if <f> and <7 are complex scalar fields 
in a theory with the symmetry breaking 17(1) x 17(1)' — ► U{ 1), the spectrum contains both 
cosmic strings and non-topological solitons. For appropriate couplings, the strings carry 
trapped charged zero modes which, when excited, act as a supercurrent. 10 In this case, 
NTS’s can be spawned in the decay of charged superconducting strings, when the charge 
condensate separates from the underlying topological vortex. 8 ’ 17 As a final speculative 
possibility, in some models non-topological solitons might form in the gravitational collapse 
of density perturbations, due to unstable spatial fluctuations in the classical a field. 

Thus, there are a variety of ways in which non-topological solitons may be born in 
the early universe. However, in the temperature range Tq > T > Iq, they are not the 
lowest free energy states of given charge, and are vulnerable to evaporation. To study 
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this process, we review the properties of NTS’s in. equilibrium at finite temperature in 
the next section. We shall confine our attention to temperatures at least slightly below 
the Ginzburg temperature Tq, so that finite temperature corrections to the structure of 
NTS’s 10 may be neglected. 


IV. Non-topological Solitons in Equilibrium 


We are interested in the question of whether NTSs formed in a primordial phase 
transition could have survived to contribute significantly to the energy density of the 
Universe. First consider NTSs in thermal equilibrium, where a distribution of lumps of 
all charges Q coexist with free <f> particles. In kinetic equilibrium, the number density of 
NTSs of charge Q at temperature T is 


n Q( T ) = 9Q 


( M(Q)T \ ,/! 

\^Q)-M(Q) 1 

( *■ ; " p 

T 


( 8 ) 


where gQ is the internal partition function for a NTS of charge Q, fi(Q) is the NTS 
chemical potential and M(Q), defined in Section II, is the NTS mass. In Eqn. (8), we 
have made the reasonable assumption that NTSs are nonrelativistic, M(Q) T, and 
nondegenerate, n(Q) C T. At temperatures T <£ m#, the distribution of free <f> particles 
obeys an expression analogous to Eqn. (8). 

Chemical equilibrium between the relative number densities of ^’s and Q’s is estab- 
lished if the accretion and evaporation reactions, (Q+ 1) «-» (Q)+<f>, occur rapidly compared 
to the expansion rate. In this case, the chemical potentials obey 


H(Q + 1) = H(Q) + nit) . (9) 

Inverting Eqn. (8) and substituting in Eqn. (9), we obtain the familiar Saha equation 

uq^ = g_Q ±1 / M{Q + 1) \ 3/2 f _ 2 tt \ 3/2 j q , T 

n Qn<t> 9Q9* V M{Q) ) \m4Tj ’ 1 
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where Iq is the NTS binding energy per charge [Eqn. (6)]. Unless the minimum stable 
NTS charge is very small, Qmin — 1> the reactions above are in fact not sufficient to 
establish chemical equilibrium. If, as an initial condition, the universe contains only <f> 
particles, accretion alone obviously cannot generate solitons. The missing ingredient is the 
fusion reaction, <f> + <j> + ... + <f> — * Qmin • 

Typically, one expects chemical equilibrium to be maintained until either the fusion 
or the accretion/evaporation rates become slower than the expansion rate. In principle, 
the proper way to study this freeze-out process is to integrate the (infinite set of) coupled 
rate equations for the number densities of all charges and free particles, i.e., numerically 
solve the Boltzmann equations for the system. We do not follow this approach. Instead, 
in Appendix A, we estimate freeze-out number densities using the “poor cosmologist’s” 
approximation: the freeze-out densities are roughly the equilibrium densities at the epoch 
when the reaction rates drop below the expansion rate. 

Such an approach, however, assumes that chemical equilibrium is initially established 
and then lost. For most NTS models, this is unjustified. In particular, when Q m in is not 
very small (as is usually the case) and/or there is a potential energy barrier between the true 
and false vacua (the NTS vacuum is a metastable minimum), then, as we show in Appendix 
A, the fusion rate is strongly suppressed and can never approach equilibrium. (Essentially, 
this is analogous to what would happen in big bang nucleosynthesis if the lightest stable 
element had a large baryon number.) In the limit that fusion can be neglected, the total 
number of NTSs is a non-increasing function of time. Since new lumps are not being 
produced, it makes sense instead to follow the evolution of individual lumps of initial 
charge Q separately. (We shall also assume that accretion of lumps by lumps is rare.) We 
develop this approach in the next Section. 
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V. Evolution of Single Lumps 

To follow the evolution of individual solitons, we need to calculate the emission and 
absorption rates of massive <p particles by NTSs of charge Q. To do this, we apply detailed 
balance arguments in the usual way. 11 First, we massage the Saha equation into a more 
convenient form. Since <p is a charged field, we assume it has g# = 2 degrees of freedom. 
In addition, up to corrections of order T/m^, we expect gQ+i/gQ — 1, and we can further 
approximate M(Q + 1 )/M(Q) = 1 + 0(p/Q). Substituting into Eqn. (10), we find 

= • (n) 
z \m^,i j 

In chemical equilibrium, the rate at which NTSs emit (p's is equal to the absorption 
rate. Since the binding energy grows with an increase in charge, absorption is energetically 
favored. Thus, if the particles are sufficiently strongly interacting, we expect a geometric 
cross section for absorption, 

* = iTUR(Q) 2 =i*u(—) Q 2r/3 , (12) 

\ m 4>J 

where the NTS radius R(Q) is given by Eqn. (2). Here, /^(< 1) is the absorption efficiency, 
which depends on coupling constants in the theory; we discuss it in more detail in Appendix 
B. The accretion rate per unit volume is simply 

■E[(Q) + (f> -►(<? + 1)] = n^nqav^ , (13) 

where v $ — (T/27^m^) 1/,2 is the mean velocity of free non-relativisitic <p particles. The 
emission rate from a lump of charge Q + 1 is proportional to the number density of lumps 
times the evaporation rate, r evap , 

R[(Q + 1) — ► (Q) + <p\ = n Q+l' r evap • (14) 
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The two rates must be equal in chemical equilibrium; using Eqns. (11-14), the evaporation 
rate from a NTS is 

= QWcxpl-Iq/T] . (15) 

7r m* 

Although the evaporation rate was derived from the assumption of chemical equilib- 
rium, Eqn. (15) applies more generally; in particular, it is valid whenever the <j> particles 
are kept in kinetic non-degenerate equilibrium [Eqn. (8)], e.g., by collisions with the am- 
bient plasma. (The power of detailed balance is that it applies to situations outside of 
chemical equilibrium.) Thus, the accretion and evaporation rate from a lump of charge Q 
is 

^=Tab.-r evap =-U0 2 —Q 2 r 13 V~2% g) -exp [-Iq/T\ . (16) 

Eqn. (16) describes the evolution of an individual soliton in a thermal bath of <j> particles. 
Before applying it, we must expose two additional hidden assumptions. 11 i) As it evaporates 
or accretes, an NTS is cooled or heated with respect to the environment; we assume, 
however, that it is kept in good thermal contact with the surrounding medium, ii) We 
assume the exterior <f> gas is dilute, so that the flow of particles into or out of the lump is 
not inhibited. 

For the remainder of this paper, we suppose there is a primordial asymmetry of <f>'s 
over anti-^’s, rj * = [n* — n*]/n 7 , where n 7 = 2 AT 3 fir 2 is the photon number density 
at temperature T. We shall further assume that <f> and are stable, but that they can 
annihilate into a pair of massless particles, say, neutrinos or photons. (These ass um ptions 
will be relaxed in paper II.) The latter assumption guarantees that annihilation takes place 
inside as well as outside solitons. For stable the Q asymmetry can be written as 

77 * = 2.5 x 10- 8 n*A 2 (GeV/m*) , (17) 
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where h is the Hubble constant in units of 100 km/sec/Mpc, and d,p is the <p energy 
density in units of the critical density. Requiring that (p particles and charges do not 
overdose the universe, d^h 2 < 1 , yields an upper bound on 77 ^ as a function of We 
define j as the fraction of (p particles in the NTS phase, so the number density of free <p 
particles is n^,{T) — 77^(1 — j)n 7 . The initial value of j depends on the efficiency with 
which NTSs are formed; e.g., in case (ii) of Section III, we expect ji ~ Pfalte < 0.5, while 
for case (i), Witten envisioned ji ~ 0.9. In general, however, j is a temperature dependent 
function which increases (decreases) during the accretion (evaporation) epoch. With these 
definitions, the expression for the evolution of a lump becomes 

§ = fi.92^(i - j) , (18) 

dl TV TT~l 0 y ^71^ J 

where f(Q) is the binding energy per <p mass [Eqn. (7)]. 

The evolution of a non-topological soliton is determined by three temperature scales. 
The evolution begins at the formation temperature Ti, which is of order the Ginzburg 
temperature Tq or smaller, depending on the formation mechanism. Once formed, the 
competition between the two terms in the brackets of Eqn. (18) determines the fate of 
the soliton. These terms are equal at the “turnaround” temperature T*, which marks the 
transition between evaporation and accretion phases. Finally, the freeze-out temperature 
Tp is defined as the temperature at which both the emission and absorption rates fall 
below the expansion rate, (\Q\/QS)t f = 1* Consider a model with ~ m^. From 
Eqns. (17) and (18), at T ~ evaporation dominates over accretion due to the small 
asymmetry term. Thus, only those charges which are large enough to survive down to T t 
or Tf are present today. There are two exceptions to this general picture. 1) If ^ is an 
unstable particle with a lifetime shorter than the age of the universe, the asymmetry is 
not constrained by Eqn. (17), and may be large. 2) If <j> is stable, but only annihilates into 
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particles more massive than f(Q)m<t>, then annihilation inside NTSs is suppressed at low 
temperature. In this case, NTSs can grow by accreting both <f >' s and <^’s, and the absorption 
rate is independent of r We discuss both of these scenarios in paper II. 

To study lump evolution, it is convenient to introduce the dimensionless temperature 
variable, 

x = m^/T . (19) 

During the radiation-dominated early universe, t — Q.ZgZ'^^mpi/rrifyx 2 (where <7. is 
the effective number of relativistic degrees of freedom), so the evolution equation can be 
rewritten, 

^ (^) . ( 20 ) 

From the discussion above, there are several possible fates that may befall a NTS with 
stable a) death by evaporation; b) evaporation followed by freeze-out, Xi < x p < x t ; c) 
evaporation followed by accretion and then freeze out, x< < x t < xp; d) accretion followed 
by freeze-out, x t < *i < *Fj and e) frozen at birth, x,- > xy. Which of these fates obtain 
depends on the mass scale, the NTS charge, the asymmetry parameter, and the formation 
temperature. 

The turnaround temperature 2*, at which accretion takes over from evaporation, is 
given implicitly by 

lnl.92 ^(l-j) = 1 . 5 lnx t -/(Q)x t . (21) 

Eqn. (21) may be derived either by equating the accretion and evaporation terms in Eqn. 
(20) or by setting tiq = tiq + i in Eqn. (11). In Figure 2 we show x t = rn^/Tt as a function 
of 77^(1 — j), for several values of f{Q). Two trends are obvious. For fixed binding energy 
f{Q), turnaround occurs ‘later’, i.e., at larger x t , as 77^ decreases, since the supply of free 
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<j> particles is dropping. Thus, from Eqn. (17), for a large mass scale m^, the evaporation 
phase is relatively long in units of the mass scale. Second, NTSs with lower binding energy 
per unit mass turn around later, as expected. For models with p < 1, f(Q) — * 1 as Q — ► oo, 
so lumps with large charge spend less time in the evaporation phase than their smaller 
brethren. (The structure in Fig. 2 at very large asymmetry, rj > 1, is only relevant if <f> is 
unstable, so we postpone discussion of it to paper II.) 

Now consider the freeze-out temperature Tp. During the radiation-dominated epoch, 
the expansion rate is FT = 1. 67 gV 2 T 2 /mpi = 1.67 gl^ 2 (rrt^/mpi)x ~ 2 . To obtain an analytic 
estimate of xf, we consider case (b) separately from (c) and (d). For case (b), freeze-out 
occurs during the evaporation phase, when the absorption term can be neglected; using 
Eqn. (20), we find 


xp = 


, (®f < *t) • (22) 


TTTvT In ( - ) + 43.9 - In (£± 7 ) - (l - ln<? 

HQ) [ \2.62 S l 12 ) 'G tVl V 3 / J 

Note that as the mass scale or the charge Q increases, freeze-out occurs earlier, i.e., 
at smaller x. Also, for large relative binding energy /(Q), NTSs freeze out at early times 
since more strongly bound objects evaporate less efficiently. 

For fixed mass scale, the freeze-out epoch xp grows as Q drops (unless 1, see 

Appendix B). Therefore, for small enough charge freeze-out will occur after turnaround. 
We denote the charge at the transition point, xp = **, by Q eq \ for Q < Q eq , we have 
xf > *t> *.e., case (c) or (d) above. From Eqns. (17), (21) and (22), 

<?., = (2.2 x (^) . (23) 

\ x t 9* / 

For p = 3/4 and 1, Q eq scales as m^ 4 and m^ 8 respectively. The dominant dependence 
in Eqn. (23) resides in the first two terms. Thus, approximately independent of p, for 
T7i0 > O(10 4 ) GeV, we have Q eq < 1. Since Q eq is the largest charge which can have an 
accretion phase, this mass scale marks an important boundary. 
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Given that NTSs are generally born in the evaporation phase (except for case (d)), we 
would like to know how large the initial charge must be in order to survive to the accretion 
or freeze-out epoch. Approximately integrating Eqn. (20) during the evaporation phase, 
the smallest charge which survives is 


Q.M - (4.6 x io»)3/(»-ap) ( 


—3/(3— 2p) ( ( l-2p/3)Ul3 2 e-f^ x< 

(1 + f(Q)xi)gl /2 


3 /( 3-2 p) 


(24) 

where we have assumed xp,x t 3> ®,- and have used a standard approximation for the 
exponential integral. In arriving at Eqn. (24), we have also made use of the fact that, for 
Q ^ 1, f(Q) is approximately independent of Q. For models with p — 3/4, the surviving 
charge scales as m^ 2 , while for solitons with p = 1, Q, scales with m^ 3 . 

From Fig. 2, we see that typically x t ~ C?(10-100). Thus, comparing Eqns. (23) and 
(24), we find Q,(x,) > Q eq if /(<?)*» < 17 + ln(m^/GeV). That is, unless NTSs form at 
a temperature well below /(Q)m0, cases (c) and (d) are eliminated: all those that would 
accrete before freeze-out in fact do not survive. If NTSs are created in a phase transition, 
we generally expect them to form at a temperature within an order of magnitude of their 
mass scale, .1 < < 10, for which we find Q , > Q tq . This can be circumvented, on 

the other hand, if NTSs form well after the phase transition is completed, for example by 
pinch-off in the asymmetric vacuum model, or by the decay of cosmic strings. In these 
cases, Xi "%> 1, Q t (xi) can be substantially reduced, and an accretion phase is possible. 
For example, with p = 3/4 and f(Q) = 1 we obtain, for Zi = 1 and se< = 10, respectively, 
<?.(!) ~ 10 33 (GeV/m*) 2 and Q.(10) ~ 10 24 (GeV/m*) 2 . 


For NTSs formed at relatively small values of ®,-, then, the only hope of survival is 
to start with sufficiently large charge, Q > Q t . However, this runs into two potential 
difficulties. First, for most of the formation mechanisms discussed in Section III, large 
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charges are exponentially suppressed and thus may be rather rare. As a result, their 
contribution to the present energy density of the universe is likely to be small. Second, 
it is difficult to imagine forming lumps with a charge larger than that contained in the 
particle horizon at the time of formation. The total charge inside the horizon is 


^ , . ..G^x 3 ( GeV\ 4 


(25) 


For example, for p = 1, with f(Q ) ~ Xj ~ 1, we have Q,( 1) > Qh( 1) for m# > 10 -4 
GeV; that is, solitons formed at early epochs (ij ~ 1) in these models are never large 
enough to survive ( Q > Q,) if their formation is limited by causality. Only if they form 
at a sufficiently late epoch, since Q t depends exponentially on f(Q)xi , can this problem 
be overcome. For models with p = 3/4, on the other hand, at Xj = 1, we find Q«(l) < 
Q//(l) if < 10 fl GeV. In this case, models with sufficiently small mass scale can 
have surviving NTSs which form early. Again, at later formation epochs, the situation 
can change substantially; e.g., for x< = 10 and p = 3/4, we find Q # (10) < (?h( 10) for 
m* < 10 12 GeV. 

In Figure 3 we summarize our results for the case in which there is a primordial charge 
asymmetry. Here we have chosen /^/? 2 = 1, g m = 100, Q^h 2 = 1, f(Q) = 1 and the two 
cases p = 3/4 and 1. As a function of mass scale, we compare the chaxge inside the 
horizon at formation, Qh(* «)> with the smallest charge that survives evaporation, 
for Xj = 1 and 10, obtained by numerically integrating the rate equation, Eqn. (20). For 
models with Q#(x,-) > Qh(xx), even horizon-sized solitons disappear by evaporation. For 
f(Q)xi = 10, both models have a window of survival if < 10 12 GeV. We also show the 
charge Q eq at which freeze-out and turnaround coincide, xj ? = x*. Recall that if Q > Q eq , 
the soliton will freeze out before turnaround; in the opposite case, there may be accretion 
before freeze-out. However, for all models and mass scales, Fig. 3 shows that Q eq << Q, 
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so that these charges die before they have a chance to grow. In Fig. 4, we show results of 
the numerical evolution of Q(x) for a sample model with p = 3/4 and m <p — 10 4 GeV, in 
the neighborhood of <?,(1) — 10 28 . For a very narrow range of Qi, the value of the charge 
at freeze-out depends sensitively on the initial charge. Outside this range, the final charge 
is either zero (for Qi < Q,) or ~ Qi (for Qi > Q t ). 

As noted above, for m $ < 10 4 GeV, it is possible to have xp > so that an 
accretion phase is possible before freeze-out occurs. Thus, even if all solitons formed at 
Xi ~ 1 evaporate away, a later period of accretion starting from individual charges can 
build up new solitons before freeze-out occurs. 13 As noted in Section IV, this is only feasible 
for very small values of the minimum charge, Q m i n ~ 1, so that free <j> particles can act as 
seeds for new lumps. We return to this possibility in Appendix A. 

VI. NTS Boiling 

Although NTS’s may be formed at a phase transition at T c , they are not the lowest 
free energy phase until the universe cools below It. In the last section we studied the 
struggle for survival of NTSs against surface evaporation. In this section we discuss the 
boiling of NTSs, i.e., bubble nucleation of the lowest free energy phase (massive <j> particles 
in the true vacuum) throughout the volume of a NTS. 

Imagine heating up a NTS from zero temperature. Surface evaporation starts at 
temperatures just above T t (or Tp, if T t < Tf). Bubble formation is negligible for temper- 
atures just above It, due to the energy cost associated with the bubble surface. However, 
at higher temperatures (T ~ Tq), bubble nucleation throughout the volume may domi- 
nate the evaporation process. The probability of bubble nucleation is very sensitive to the 
surface tension of the bubble wall, <Tb , which is determined by the mass scale and coupling 
constants. NTSs can survive the boiling phase if <Tb is greater than a critical value derived 
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below. 


Following Ref. 12, we calculate the thermodynamic work expended to create a bubble 
of free <f>’s inside a NTS: 

W = — —r 3 (Pb ~ Pnts) + 47r cryr 2 , (26) 

w 

where r is the radius of the bubble, Pb is the pressure inside the bubble, and Pnts is the 
pressure inside the NTS. W is maximized at the critical radius r c = 2 <rb/(Pb — Pnts )• 
Bubbles with radius r < r c shrink away while those with r > r c grow. The nucleation 
rate, T, is determined by the abundance of critical size bubbles, 

r ~ T A e~ Wc/T , (27) 


where 


W c = 


16t r 


(28) 


3 ( Pb-PNTS ) 2 

The soliton is assumed to be in pressure equilibrium with its surroundings. Therefore, 
Pnts Is just equal to the exterior thermal radiation pressure. The pressure inside the 
bubble, Pb, has the same thermal component, plus the pressure due to the massive <f> gas, 
P^. Hence, 

Pb - Pnts = P* = (2/T 3 )^ 2 m 3 J 2 T 5 / 2 e xp(-I Q /T) , (29) 


and 

Wc _ 8 tt 4 <rj 2lo jT 

T 3 mjT 8 

At temperature T, the characteristic number density of bubbles, n&, is of order Tt, 
where t is the expansion timescale, so ~ 0.2g^ 1 ^ 2 mpiT 2 exp(—W c /T) . If rib » w NTS > 
the charge density inside the soliton, the NTS phase is so far out of thermodynamic 
equilibrium that the probability of remaining in such a phase is negligible. (The NTS would 
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likely not have been formed, let alone survived.) The condition rib < tints translates into 
a constraint on the surface tension crj, > where 

1 1/3 




(31) 


Here we have assumed a spherical NTS lump of charge Q and radius Rq given by Eqn. 


( 2 ). 


On the other hand, if <7& is too large, boiling may never occur. Boiling will only 
dominate surface evaporation if the surface area of the NTS lump is smaller than the 
surface area of all the bubbles inside the lump. Boiling is more effective than surface 
evaporation if crj, < <r moz , where <r mit is given implicitly by: 

<rTnax{x,Q,m<i,) = .57m^ 


3 exp(-2f(Q)x/3) 


x 

(32) 

Both <T m i n and <j m ax depend on the specific NTS model. In Fig. 5, we plot <7 mtn (x) 
(lower curve) and <r maz (x) (upper curve) for different values of m^, where we set Q = 
Qhor{tti^,x = 1) [Eqn. (25)]; we have chosen f(Q) = 1, except for the case = 1 GeV, 
where we set f(Q) = .05, characteristic of strange matter. We have also set p = — 1 , 

since the dependence on /? (~ 0(1)) is only logarithmic, and the difference between p = 1. 
and p = .75 is unnoticeable. For < <r m i n , there is copious boiling and NTSs disintegrate. 
Between the two curves, <r m in < &b < <r max , boiling dominates over surface evaporation, 
and for <rj, > tr max boiling can be neglected. 

From Fig. 5, we see that, for models with f(Q) ^ 1, boiling is important if x,- ~ 1. 
In these models, the surface tension crj, must be larger than about .3m^ to avoid boiling. 
Alternatively, for smaller <76, NTSs can survive if they form later. More generally, the 
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efficiency of boiling is large for f(Q)xi < 1. For example, for strange matter, f(Q) — 0.05, 
X{ — 10, and strange lumps can only survive if 12 trj, > .006m^ (mucli larger than the 
theoretically expected value 4 ). 

VIII. Conclusion 

We have discussed the cosmic evolution of non-topological solitons, developing a gen- 
eral analysis which applies to virtually all classes of models, and thus separating general 
features from model-specific results. We have found that, for a wide range of parameters 
and models, primordially formed NTSs of moderate charge disappear quickly, either via 
surface evaporation or bubble nucleation, leaving behind a distribution of free massive par- 
ticles. For a range of mass scales, depending on the scaling parameter p, there is a window 
of survival between the smallest charges which survive evaporation, Q, , and the largest 
that can be causally formed, Qh • This window, in turn, depends on the epoch of NTS 
formation; in scenarios with delayed NTS formation, the chances of survival are enhanced, 
and, if < 10 4 GeV, sufficiently small charges may then undergo an accretion phase. 
Our results pertain to models constrained to have small charge asymmetry by the existence 
of stable <(> particles. If this constraint is relaxed, the results can change considerably (cf. 
paper II). 
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Appendix A 

In this Appendix, we consider a system of interacting <f> particles and solitons which 
reaches chemical equilibrium for a time in the early universe. Recall that this might occur 
in two ways. 1) If Q m in ~ 1, chemical equilibrium can be established by the accretion 
and evaporation reactions, (Q + 1) «-» (Q) -j- (f>. (This possibility was recently studied 
by Griest and Kolb 13 in the context of a particular toy NTS model with p = 3/4.) 2) 
In most NTS models, unless the coupling constants are specially chosen, Q m *n 1. In 
this case, the fusion and dissociation reactions ^ + <f> + ... + <f> <-» (Q) are also required to 
establish complete chemical equilibrium. However, at the temperatures at which fusion 
can be effective in generating new solitons, T < T t , these reactions are likely to be strongly 
suppressed compared to the expansion rate, for reasons given at the end of this Section. 

In both cases, in chemical equilibrium the chemical potentials are related by 

H{Q) = Qfi(<f>) , (A.l) 

where axe given implicitly by Eqn. (8). We define the mass defect by 

B q = Qm* - M(Q) = - aQ '"») = QmJ(Q) , (A. 2) 
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and note that f(Q) -+ f(Q) [Eqn. 7] as p -► 1 or as Q -* oo, independent of p. Then, from 
Eqns. (1), (8), (A.l) and (A.2), the number density of solitons of charge Q in chemical 
equilibrium can be written 

_ / s3(Q-1)/2 

n “ iT) = (^f) eBq ' T ' (A3) 

It is convenient to introduce the charge fraction of solitons of charge Q, Yq — uqQ/N , 
where the total charge density is given by 

OO 

N = ^ Q n Q — r f<t > n 7 ■ (^- 4 ) 

Q min. 

Note that, in general, the charge fraction is different from the mass fraction. From Eqns. 
(A.3) and (A.4), we have 

Yo = V /J l (A. 5) 

J 

and the constraint 

E Y <,-i ■ u- 6 ) 

Qtnin 

The charge fraction is obtained by simultaneously solving Eqns. (A.5) and (A.6). 

The qualitative behavior of the NTS population is relatively simple to understand. 
At high temperature, the relative abundance of NTSs with large Q is suppressed by the 
small asymmetry factor ~ p^ -1 . At lower temperature, the exponential and Q 1+ ( 3 p/ 2 ) 
terms favor the formation of large Q objects. The trend from small to large Q shifts at the 
‘turnaround’ temperature defined in Section V: in the Saha equation, we find nq = tiq+i 
at T = Tf. Since f{Q) is approximately constant for large Q, x t is roughly (^-independent. 
Thus the turnaround is rapid: slightly above T t , large charges are strongly suppressed, 
while slightly below it, they are strongly favored. As the temperature drops through T t , 
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large charges quickly build up until the available ‘fuel’ of free particles is depleted, thereby 
shutting the system out of chemical equilibrium. 

At high temperature, x < x t , Yq falls rapidly to zero with increasing Q. Consequently, 
in this regime, only the first few terms contribute significantly to the sum in Eqn. (A. 6), 
and the equilibrium charge fraction may be easily calculated. Thus, if x F < x t , i.e., 
if freeze-out occurs before turnaround, the final abundance Yq[x f (Q )] can be estimated 
analytically. As discussed in Section V, this is always the case for m $ > 10 4 GeV. At 
x < xt, we expect the total NTS charge fraction to be small, j <GL 1. To zeroth order of 
approximation, we can then estimate Yq by setting j = 0 in Eqn. (A.5). 

Alternatively, we can obtain a more transparent expression for Yq(x f ) without ap- 
proximation. Setting g<f> = 2 as in most cases of interest, the freeze-out abundance for 
x F < Xt may be written 

Y Q (x F ) = 10 7 (a*F) 1/2 Je ) (n,ft 2 r 1 <? 1+(2 >’' 2 > (4.8 x 10 "*)« 


(1 -j) ( 4.4 X 10 


hQ)/f{Q ) 


3/2 


s : /2 


/GeV\' t(te)rt( ® 1 

Kw) 


Q 


(A.7) 


where x F (Q) is given by Eqn. (22), f(Q) by Eqn. (7) and f(Q) by Eqn. (A.2). If we 
approximate f(Q) = f(Q), valid at large Q, this may be rewritten using Eqn. (23), 


Yq{x f ) = 10 7 


(az F ) 3 / 2 gQ 


gl+(3p/2) 



(A.8) 


This gives the freeze-out abundance for Q > Q eq and thus x F < x t . The case x F ~ x t 
must in general be handled numerically. 

We finish this section by discussing the fusion of free <f >* s into a soliton for models in 
which Qmin is not small. This rate must be faster than the expansion rate in order to 
maintain chemical equilibrium. We shall show that this is unlikely at temperatures below 
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m <p. We estimate the rate by assuming that a soliton is formed when Q > Q m in <f> particles 
get together in a fluctuation with a charge density comparable to the NTS’s charge density. 
This criterion can be made rigorous by considering the effective c potential U(cr,q) in a 
background of non-zero charge density q. At a charge density q ~ qNTS, the unconfined 
phase a = (To becomes degenerate with the NTS phase and at higher densities the NTS 
phase is preferred. If the vacuum potential U(<t, 0) has no barrier (see Fig. 1), then a 
soliton forms whenever a sufficiently dense charge fluctuation occurs. However, if there 
is a barrier between the NTS state and the unconfined phase (dashed curve in Fig. 1), a 
charge fluctuation is not sufficient to form an NTS: in addition, the cr field must tunnel 
through the barrier, leading to an extra suppression of the fusion rate. We will ignore this 
factor, so the estimate below will yield an upper bound on the formation rate. 

For simplicity we focus on the case p — 1; the results may be easily generalized to 
other models. Consider a fluctuation region with fixed volume V’ fmn = (47r/3)/3 3 m^ 3 Q min , 
the size of a NTS with charge Q m in, containing an average charge Q = qV m i n , where the 
mean charge density q ~ (2Arj < f > /7r 2 )T 3 . The probability of finding a charge Q in a region 
■with average charge Q is of order 


P{Q,Q) - (2ttC?) 1/2 exp 


( q-q r 

2 Q 


dQ . 


(A.9) 


Here we have assumed Q is of order a few or more, so that the Gaussian approximation 
to the Poisson distribution is reasonably accurate, and have ignored antiparticles. We 
are interested in a fluctuation with charge Q > Qmin- At the temperatures of interest, 
T < m^, we can safely assume < (m^//?T) 3 , so that Q m %n > Q • Thus, the probability 
of a fluctuation with charge large enough to make a soliton is 

P(Q > «™„) a (^) • M-10) 
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The number density of minimum charge NTSs is n to i ~ P(Q > Q m in)/Vmin , and their 
abundance relative to free <j) particles is approximately 


n BOl ^ 

( * 3 V' 3 ... 

Qmin 

( * 3 'll 

n 4, 

\P*Q\ninV4>J 

[ 2 



Thus, if Q m in is large and/or the asymmetry 77 ^ is small, the formation of solitons from 
thermal fluctuations is strongly suppressed at temperatures below m $ [x > 1). 

Appendix B 

The factor measures the probability of absorption into a NTS of an* incoming 
particle (f>, relative to its geometric cross-section. Since the mass term for <f> generally 
drops abruptly (and monotonically) as the NTS surface is approached from the exterior, 
an incoming <f> particle effectively sees a square well potential inside the NTS. Since there is 
no potential barrier, <j> particles have no difficulty entering the NTS, but they can as easily 
escape at the other side. The incoming <f >' s (massive outside the NTS) have energies higher 
than the binding energy of the massless <f> gas inside the NTS. Therefore, is determined 
by the probability of an incoming 4> being trapped by scattering down to lower energies. 

The probability of at least one collision inside a soliton of radius Rq is given by 
Pc — 1 — exp(— .Rq/A^) , where is the mean free path of </>'s inside the NTS. Thus, 
= (n,<T^) -1 , where n, is the number density of scatterers and is the cross section for 
^-scattering. If Rq A^, the incident <f> will have enough scatterings to be trapped, and 
ft = p e ~ 1 . On the other hand, if Rq < A^, then p e ~ Rq/ A^ < 1 . In this case, there 
will typically be at most one scattering, and fa is given by the probability of scattering 
and losing enough energy in one collision to be trapped. 

If a <f> particle has initial kinetic energy E outside the NTS, it will have energy E{ = 
E + mj, inside the NTS. (Here, we are assuming <f> is massless inside an NTS.) The particle 
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is captured if it scatters down to a final energy inside the NTS, Ef < m In the case 
of multiple scattering, capture generally occurs, since the incident energy E ~ T < m^. 
In the single scattering regime, X# < Rq , the capture probability is model dependent. 
If we consider isotropic cross sections, on average half of the incident energy is lost per 
collision. If the target particles responsible for capture are not <f> particles, the condition 
for trapping in one scattering is roughly E < m^. However, if the dominant scatterers are 
<f>' s, the incident energy must be lower; otherwise, the initially bound <f> will be scattered 
out of the NTS, with no net gain in charge. In this case, the condition for single scatter 
capture is approximately E < Jq, and thus T < Iq. 

The mean free path A^> is also model dependent. In the simplest NTS models, the 
possible interactions are either or <f>-<r scattering. In some models, although the cross- 
section for <f>-<r scattering is appreciable, the density of tr particles inside the NTS is small 
(e.g., if a is very heavy). In these models, scattering can occur either through a <f> A 
term in the Lagrangian, or through the exchange of a c particle. 

For example, consider a model where <f>-<f> scattering through a g<f>* interaction dom- 
inates, with cross-section a, = g 2 1 Z2irE}. If the <j)' s inside the NTS are non- degenerate, 
the density of scatterers is n, = 3Q/(4ttRq). Then 


ft = 1 — exp 




1— 2p/3' 


(B.l) 


where 


_ ( 128w 2 (3 2 (l + *n 
Qm V 3 g 2 x 3 ) 


3x 3/(3-2 p) 


(B. 2 ) 


In this case, ~ 1 for Q » Q m (multiple scattering), and ~ ( Q/Q*y~ 2p ^ ^ or Q < 
Q m . At x = 1, for models with p = 3/4, we find Q 9 ~ 10 *((3/g)*, while for p = 1, 
Q , ~ 10 9 (/3/g) 9 . Thus, the approximation /^ = 1 is valid for large charges, but for small 
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charges the scattering probability is suppressed. Consequently, small charges freeze out 
earlier; for Q < Q„, 


xf — x F (U = 1) - ^1 - In 


(B.3) 


and 


Qeq — (q ^ Qtqifj — l) > 

(BA) 

Q. = (£) Q.(U = l) • 

(B.5) 


In principle, the distribution of NTS’s with charge Q Q* may be frozen out from the 
beginning (xp < Xi) and thus survive. For these objects, the mean free path is much 
larger than the NTS size, so the timescale to repopulate the evaporated tail of the interior 
<f> distribution is longer than the expansion timescale. In this case, the initial distribution 
of small NTSs is preserved, and may contribute significantly to the energy density of the 
universe. On the other hand, if <f> is kept in thermal equilibrium by scattering with other 
particles, the mean free path is much shorter and we expect ~ 1 down to smaller values 
of Q. 
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Figure Captions 

1. Schematic a (confining field) potentials for theories with NTS solutions, with a pre- 
sumed NTS phase at <r — 0 and true vacuum <ro = 2.5 (arbitrary units). For the 
dashed curve, the ‘NTS vacuum’ is metastable in (the absence of charge), while for 
the solid curve it is unstable. 

2. The turnaround epoch x* = m^/T as a function of the asymmetry parameter r], for 
different values of the binding energy per unit mass, f(Q) = 0.05, .25, .36, .5, and 1. 

3. The minimum charge which survives evaporation, Q,(xi), is shown as a function of 
the mass scale for the two cases p = .75, 1, and for ® t - = 1 and 10. Also shown are 
the charge inside the horizon, Qh(*») and the charge below which an accretion phase 
is possible before freeze-out, Q eq . 

4. A sample numerical integration of Q(x), for the parameters f(Q) = 1, = 10 4 GeV 

and p = 3/4. In this example, Q, ~ 10 28 . 

5. The surface tension <r m i„ and cr mai as a function of x, for mass scales = 
10 2 ,10 5 ,10 15 GeV, for f(Q) = 1. Also shown is the strange matter case, f(Q) — 0.05, 
m<£ = m n ~ 1 GeV. Boiling dominates evaporation for cry, < <r max . 
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